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We investigate theoretically the eigenmodes and the stability of one and two arbitrary fractional vortices 
pinned at one and two K-phase discontinuities in a long Josephson junction. In the particular case of a single K- 
discontinuity, a vortex is spontaneously created and pinned at the boundary between the and K-regions. In this 
work we show that only two of four possible vortices are stable. A single vortex has an oscillatory eigenmode 
with a frequency within the plasma gap. We calculate this eigenfrequency as a function of the fractional flux 
carried by a vortex. 

For the case of two vortices, pinned at two K-discontinuities situated at some distance a from each other, 
splitting of the eigenfrequencies occur. We calculate this splitting numerically as a function of a for different 
possible ground states. We also discuss the presence of a critical distance below which two antiferromagnetically 
ordered vortices form a strongly coupled "vortex molecule" that behaves as a single object and has only one 
eigenmode. 

PACS numbers: 74.50.+r, 85.25.Cp 74.20.Rp 
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I. INTRODUCTION 

Vortices in long Josephson junctions (LJJs) usually carry 
a single magnetic flux quantum <t>o and therefore are often 
called fluxons. The study of dynamics of fluxons has been 
attracting a lot of attention during the last 40 years because 
of their interesting nonlinear natur e) 1 : 2 ! 3 as well as because of 
potential applications^. 

Recently it was shown that one can fabricate so-called 0- 
7t-LJJs, i.e., LJJs consisting of alternating regions with posi- 
tive (0-part) and negative (ft-part) critical currents. Such junc- 
tions can be fabricated using superconductors with anisotropic 
order parameter which changes sign depending on the direc- 
tion in fc-space (e.g. d-wave order parameter symmetry)^ or 
with an oscillating order parameter (e.g. with a ferromagnetic 
barrier)&2J2. It was shown that at the boundary between a 
and a K part a new type of vortex carrying only half of the 
flux quantum may existii*i2Ji Such vortices (naturally called 
semifluxons — ) were observed experimentally in several types 
of 0-7T,-LJj9i£ii£iIii&. They may appear spontaneously and 
correspond to the ground state of the systemi2i22iSl. 

The dynamics of the Josephson phase in a LJJ is described 
by a perturbed sine-Gordon equation. For the case of the 0-71- 
LJJ this equation is slightly modified^ and includes the func- 
tion B(x) which is equal to zero all along the 0-parts and is 
equal to 7C all along the ft-parts. The function Q(x) as well as 
the Josephson phase §(x,t), which is a solution of the sine- 
Gordon equation, is a 7C discontinuous function of x at the 
0-7C-boundaries. We sometimes call the 0-7t-boundaries dis- 
continuity points. 

Recently, a LJJ geometry which allows to create arbitrary 
discontinuities was suggested and successfully tested 2 ^. In 
this LJJ a pair of closely situated current injectors creates an 
arbitrary K-discontinuity (not only K = ±Jt) of the Joseph- 



son phase, with K being proportional to the current passing 
through the injectors 2 ^ 3 .. Since the Josephson phase is de- 
fined modulo 2k without loosing generality below we consider 
only < K < 2%. 

Similar to the case of a ^-discontinuity, the presence of an 
arbitrary K-discontinuity results in formation of a fractional 
vortex pinned at it 37 . In the following, the topological charge 
of a single vortex is defined to be a conserved quantity un- 
der the boundary condition imposed on the phase fj(x), and it 
is equal to the difference between the phase at x = +°° and 
x = — oo. We use the word K-vortex, to denote a vortex with 
the topological charge equal to K. Given a discontinuity — K, 
the topological charge of the vortex should be such that its 
sum with the discontinuity is equal to 2nn (n is an integer), 
otherwise the energy of the system diverg esA In fact one can 
construct not an infinite number but only four such vortices 
corresponding to « = —2, —1,0, 1, namely a (k — 47t)-vortex, 
a (k — 27t)-vortex, a K-vortex, and a (k + 27t)-vorte54 2 4 

It is shown below that only two of these vortices are 
stable: a so-called direct +K-vortex and a complementary 
(k— sgn(K)27t)-vortex. By definition the complementary of 
a complementary vortex, gives again a direct +K-vortex. The 
only exception is the case K = 27t«, for which there exist three 
stable solutions: a +27t-vortex (fluxon), the flat phase state 
(zero phase), and a — 27t-vortex (antifluxon). The complemen- 
tary "vortex" for both ±27t-vortices is a constant phase state 
and, vice versa, for the constant phase state the complimen- 
tary vortex is either a fluxon or an antifluxon. In the majority 
of situations this can be distinguished due to the conservation 
of the topological charge. 

Since, for K ^ 2%n, K-vortices are pinned, they might have 
an eigenmode corresponding to the oscillation of the flux 
around the discontinuity point. Imagine that we apply a small 
uniform bias current through a LJJ containing a fractional vor- 
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tex. The current creates a driving (Lorentz) force which bends 
the vortex. The driving force is compensated by the elastic (re- 
turn) force which arises due to vortex deformation. If we now 
suddenly remove the bias current, the elastic force makes the 
vortex oscillate around the equilibrium shape with the eigen- 
frequency coo, provided the damping is low enough (like in 
artificial O-K-LJJs^ 2 .). 

We note that an integer fluxon does not have any non-trivial 
eigenfrequenciesS. Once deformed or perturbed, a fluxon 
recovers its shape exponentially without any internal oscil- 
lations. When the bias current is applied, the fluxon starts 
moving in a certain direction (which depends on the polarity 
of the fluxon and the bias current) which formally corresponds 
to coo = (dc motion). 

In this paper we investigate the stability and the eigen- 
modes of the vortex states formed in a LJJ with one or two K- 
discontinuities. The paper is organized as follows. In section 
||l]we present an analytical derivation of the eigenfrequency 
of a single K-vortex and compare the obtained results with 
the results of direct numerical simulations of the sine-Gordon 
equation. Then we study numerically two coupled K-vortices 
in section ITTT1 The eigenfrequency in some limiting cases is 
derived analytically . Section lTvl concludes this work. 



II. EIGENFREQUENCY OF A SINGLE VORTEX 

For our study, it is more convenient to use the continuous 
Josephson phase fj(x,t) defined as fi(x,t) = §(x,t) — &(x)*^ A 
perturbed sine-Gordon equation which describes the dynam- 
ics of the Josephson phase in a 0-K-LJJ can be written in terms 
of n(x,t) asii 



with 



Hxx - Vtt - sin[u + Q(x)] = a/j, - y(x), 



with 



Q(x) 



0, x<0, 
-k, x > 



(1) 



(2) 



for the case of a single discontinuity at x = 0. Without loss 
of generality, we assume < K < 2k. The Josephson phase 
fi(x,t) is a continuous function of the coordinate x and time 
t, which are normalized to the Josephson penetration depth 
Xj and to the inverse plasma frequency (Op 1 , respectively; 

a = 1 / y/$~ c is the dimensionless damping parameter, p c is the 
McCumber-Stewart parameter, and y(x) = j(x)/j c is the bias 
current density j(x) normalized to the critical current density 
j c . The subscripts x and t denote corresponding partial deriva- 
tives with respect to x and t , accordingly. 

It is natural to have the boundary conditions at x = given 

by 



(3) 



The zero bias static solution of Eq. Q corresponding to a 
direct K-vortex is given by 



f(x) = 



f£(x)=f(x),x<0, 
H«(x)=K-f{-x),x>0, 



(4) 



f(x) =4 tan e 



l J*+*o) 



The constant xq still has to be determined. 

Imposing the boundary conditions Q to /J K (x), we end up 
with the expression for xq: 



xq = In tan - 



(5) 



Next we will calculate the eigenfrequency of /A The 
calculation we present in this article follows the one from 
Ref. |23 First we linearize Eq. ([0 about the solution /j k . We 
write n(x,t) = /j k + u(x,t) and substitute the spectral ansatz 
u = e )u v(x) into Eq. Q. Retaining the terms linear in u gives 
the following eigenvalue problem 



[A 2 + al + cos(v K + 8)] v = 0. 



(6) 



Note that y is assumed to be zero here. 
If we define 



, x < 0, 



v+,x>0 



then at x = the boundary conditions are 

v-(0) = v+(0), 
3,v~(0) = 9*v + (0). 



(7a) 
(7b) 



From Eq. ijfji, we know that v satisfies the equation 

vt - + ocA, + cos(/(=fx))) v ± = 0. 

The eigenfrequency corresponds to a solution v that tends to 
as x — > ±°°. According to Refs.[2(|and0fhe above eigen- 
value problem, if considered on the whole real line, has two 
independent solutions. We will use these known solutions to 
construct the eigenfunction and find the corresponding eigen- 
value for our problem, i.e., including the discontinuity. The 
solution that corresponds to the eigenfunction is given by 



where 



- (V) = e M*+*o) (tanhfx + x ) - A); 
-(x)=v-(-x), 

A= y/X 2 + aX+l. 



(8) 



(9) 



The condition that v — > as x — * ±°° requires that the real 
part of lambda, Re (A), is positive. The eigenvalue X is calcu- 
lated from A by determining that v has to satisfy Eqs. 0. 
From Eq. ilbl we obtain that 

v+(0) -v~(0) = 2 e Av (A(tanhy- A) +sech 2 y)| x=0 = 0, 
y = ±x + xq, 

from which we obtain 



A = — ( tanhxn ± \l tanh 2 xo +4 



l-sech Xq 
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FIG. 1 : Plot of the normalized eigenfrequency coo of a single K- vortex 
versus K given by Eq. lilt [continuous line] and its approxima- 
tions for K — * 0: K)o(k) w 1 — K 4 /512 [dashed line] and for K — * 2k: 
COo(k) rj j\/2jt — K [dotted line]. Symbols show the values of eigen- 
frequency obtained by direct numerical simulation of Eq. Q. 



When a = 0, the eigenfrequency coo is given by 



X = too = ±-y 4 — ^tanhxo + y tanh 2 xo + 4sech 2 xo 

(10) 

with xq given by Eq. l|5}. This can be further simplified to give 



coo(k) 



\ 




cos ■ 



4-3 cos 2 



(11) 



The plot of this eigenfrequency as a function of K is shown 
in Fig.^ For the particular case of a semifluxon the eigenfre- 
quency is goo(ji) = \ V 1 + V5 w 0.899. Note that coo(k) = 
cOo(— k) because the K- vortex and the — K- vortex have the 
same eigenfrequencies. 

The above result is the only eigenvalue of a direct fractional 
vortex. For a rigorous proof, we refer to Ref. [Iff Below we 
will show that there is no eigenfunction with at least one zero 
at finite x. 

According to the Sturm-Liouville theorem, if Eq. has 
several eigenvalues E = X 2 + aX = A 2 — 1 given by, e.g., 
E\,E2, ...,£„, then we can arrange them such that E\ < E2 < 
■ ■ ■ < E„. The eigenfunction corresponding to E„ will have 
(n — 1) zeros. 

Next, based on the first line of Eq. (|SJl, we construct an 
anti-symmetric solution with v + (x) = — v~(— x). With this 
property of v, v = at x = and Eq. Q can be fulfilled if 
and only if 



v(0)=0 or A = tanh(x ) = -cos^-J < 0. (12) 

We see that there is no solution A with positive real part, 
which is necessary to let v — > as x — > ±°°. Only at K = 2% 
we have an eigenvalue attached to the edge of the continuous 
spectrum and the corresponding eigenfunction is bounded, but 



not exponentially decaying. 2 '' At this value of K, the solitary 
wave we consider is nothing else but an integer fiuxon. Hence, 
a stable fiuxon has no internal mode. 

We have also checked our analytical results numerically us- 
ing StkJJ softwara^. We simulated a rather long L = 20Xj 
Josephson junction with a K-discontinuity at its center x = 0. 
To obtain the eigenfrequency, we applied a small bias current 
which pulls the K-vortex in a certain direction and waited to 
arrive to a stable stationary state. Then we abruptly decreased 
the bias current back to zero and observed the /j, (?) °= V(t) at 
x = 1, where V(f) is the normalized voltage. The plot V(t) 
exhibits decaying oscillations. Making a Fourier transform, 
we find a single dominant frequency which is considered as 
an eigenfrequency of the K-vortex. We repeated this proce- 
dure for different positions x and different small initial bias 
currents, but the results did not change within 1 % of accuracy. 
Repeating this numerical experiment for several K in the range 
from to 2k we obtained the values of eigenfrequencies which 
are shown as filled symbols in Fig. ^ demonstrating perfect 
agreement between the analytics and the direct simulations. 
The eigenfrequency of the complimentary (k — 27c)-vortex can 
be obtained by using 271 — K instead of K in Eq. il It . 

Note that Eq. dlOt is also valid for (k + 27Cn)-vortices 
with the topological charge larger than 2k. However, when 
|k + 27C«| > 4jc the static solution cannot be constructed 24 . 
Therefore we consider here only the cases n = 1 and n = —2. 
The zero bias static solution of Eq. Q corresponding to a frac- 
tional (2ku + k) -vortex is given by 



V K {x)=f(x),x<0, 

fJl(x) = 27tn + K- /(-x),x > 0, 



(13) 



with 



f(x) =4 tan l e 



where x„ is determined by matching the boundary conditions 
as described above. Thus, 



x n = In tan 



(Kn k \ 



(14) 



The eigenfrequency of this state is still given by dlOt with xq 
being substituted by x n . One can check by direct substitution 
of into (O that Re(X) > for « = +1 or n = -2, i.e., the 
(k+2tc) and (k — 4tc) vortices are unstable and, therefore, not 
observable. A damping cannot stabilize these "heavy" vor- 
tices. 

If one thinks about the dispersion relation co(£) = \/l +k 2 
of a LJJ, the dispersion relation has a gap (plasma gap) from 
co = to co = 1. In the presence of a vortex, there is an 
additional discrete level (eigenfrequency) situated within the 
plasma gap and corresponding to k = (localized state). This 
level is somewhat similar to impurity levels in semiconduc- 
tors. The position of this "impurity" level can be tuned elec- 
tronically by changing K. 



III. TWO COUPLED VORTICES 

Now consider two K-discontinuities (0 < K < 2tc) in an in- 
finitely long LJJ at the distance a from each other (at x = 
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Symbol Discontinuity Topological charge Name 



1 


— K 


+K 


direct 


I 


+K 


— K 


direct 


I 


+K 


271 — K 


complementary 


I 


— K 


K— 271 


complementary 


I 


±71 


+71 


semifluxon 


i 


±71 


— 7t 


antisemifluxon 


t 





+27t 


fluxon 


4 





-271 


antifluxon 



TABLE I: Notations for different types of vortices. 



±a/2). If both discontinuities have the same sign of K (e.g. 
+k), there are two possible irreducible vortex configurations: 
the symmetric ferromagnetic (FM) state JJ = (— k, — k), and 
the asymmetric antiferromagnetic (AFM) state J \= (— K, 2n — 
K)2i. If the discontinuities have different sign, e.g., +K and 
— K, there are two other irreducible vortex states: the asym- 
metric FM state |"| = (2k — k, k), and the symmetric AFM state 
J ] = (~K,+K)Sm. The details about ground states are presented 
in Ref. [24} Below we consider eigenmodes of these ground 
states. 

Note that in the notations such as 1 J \ [ the direction of the 
arrows shows the polarity of the vortex (up or down), while 
the harpoon on the left or on the right side indicates whether 
the vortex is direct or complementary. The notations are sum- 
marized in Tab. H] 

In a system of two weakly interacting vortices the eigen- 
frequency ©o splits into two different frequencies: co + , cor- 
responding to the in-phase oscillations of both fractional vor- 
tices, and co_, corresponding to the anti -phase oscillations. By 
definition adopted here the "in-phase" or "anti-phase" means 
that magnetic fields (rather than Josephson phases) of two vor- 
tices oscillate in-phase or anti-phase. According to the Sturm- 
Liouville theorem the eigenfunctions v(x) (oscillations of the 
Josephson phase) corresponding to the first and the second 
lowest eigenfrequencies have no zeros and one zero, respec- 
tively. In turn, the oscillations of magnetic field V x (x) have 1 
or 2 zeros. If one draws such V x (x) deviations qualitatively 
on the top of the magnetic field profiles ^ x {x), corresponding 
to various ground states, one can conclude which eigenfre- 
quency corresponds to in-phase and which to anti-phase os- 
cillations. 

In general, it is quite difficult to calculate the splitting of 
the eigenfrequency analytically 39 , so below we calculate the 
frequencies co± (k) for different values of a numerically. 

To obtain the data presented below we have used a tech- 
nique similar to the one mentioned in Sec. |n] To excite in- 
phase oscillations we applied a small uniform bias current. To 
excite anti-phase oscillations the bias current y(x) was applied 
so that y(x) > for x > and y(x) < for x < 0. We used 
several different y(x), e.g., y(x) «iot y(x) °= sgn(x), and the 
results were generally the same. 

Unfortunately this technique is very slow, requires a lot of 
manual work and is not very precise. Therefore, we imple- 
mented another technique for calculating the eigenvalues in 



StkJJ**. When the numerical solution decays to a stable sta- 
tionary one, we analyze its stability, by introducing an arbi- 
trary perturbation 8 to the found solution /jq, i.e., we substitute 
fi = /no + 8 into and obtain a PDE for 8. Then we discretize 
the obtained PDE and arrive to a system of N coupled sec- 
ond order ODEs. Typically for a LJJ of length L = 40Xj, we 
use the spatial discretization step Ax = Q.Q2Xj which gives us 
N = L/Ax = 2000. This system can, in turn, be reduced to a 
system of 2N ODEs of the first order, which can be written in 
a matrix form as 

8=A8, (15) 

where 8 is a vector with 2N elements and A is a constant 
matrix constructed using /jq. The eigenvalues of this matrix 
can be found using standard routines 2 - 9 . Among all 2N com- 
plex eigenvalues we chose two eigenvalues with the smallest 
imaginary part (eigenfrequency). The other eigenvalues had 
| Im(A.)| > 1, i.e., belonged to the plasma band. The real parts 
of all eigenvalues were negative (and w —a/2), indicating the 
stability of the solution. In this way, we were able to cal- 
culate eigenfrequencies numerically, with high accuracy for 
hundreds of K values in automatic fashion. The results pro- 
duced by both methods were compared in a number of se- 
lected points and were essentially the same. 

A. Symmetric AFM state 

First, let us consider a LJJ with two (+K, — k) discontinu- 
ities and with (— K, +k) vortices J } pinned at them. Since this 
symmetric AFM state is the most natural state of the system, 
in simulations we were simply sweeping K from to 2ji start- 
ing from /j = solution at K = and the symmetric AFM state 
was formed automatically. The numerically obtained eigen- 
frequencies co± (k) for different distances a between vortices 
are shown in Fig. [2] 

It was shown2i that there is a critical value of discontinuity 
Ke (a) (7t < k1 < 2ti) at which the ground state 1 J switches to 
the state [ f (k — 2k and 2k — K vortices pinned at — K and +K 
discontinuities). The inverse function for TQ (a) is denoted as 

When a > a c (k) = a c = k/2, see Fig.|2ja), the splitting of 
the eigenfrequency is clearly visible with co_ < coo < <*>+• The 
eigenfrequencies co± are calculated for < K < (a). At the 
critical value k\} (a), the state \ \ becomes unstable. As a sign 
of this instability one can see that <x>_ sharply approaches as 

K — * K c . 

The eigenfrequency of the complementary state CO^(k) is 
given by 

K>±- (k) = C0± (271 — K), K > 27C - Kc. 

Note, that in the interval of K from 27t — K e J to K C J both 1 J and 
f [ states are stable, but have different energy^. Hence, we 
have a bistability phenomenon. 

When a < a c , see Fig. |2jb), the transition from the state 
1 J to the state [ \ is smooth without hysteresis and happens at 
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FIG. 3: Eigenfrequencies C0_(k) (solid lines) and co+(k) (dashed 
lines) for the symmetric FM state \ \ calculated numerically for dif- 
ferent distances a between the vortices. The dotted lines show the 
critical values of Kc for given a at which the state \ \ ceases to ex- 
ist. Dashed-dotted lines show eigenfrequencies for a single direct or 
complimentary vortex in a LJJ with one discontinuity. 



B. Symmetric FM state 



FIG. 2: Eigenfrequencies co_ (k) (solid line) and co+ (k) (dashed line) 
for the symmetric AFM state calculated numerically for different 
distances a between the vortices in the limit of (a) weak (a > a c ) 
and (b) strong (a < a c ) coupling. For the case a > a c only the state 
1 J is shown. For a < a c both states 1 J (k < 7t) and [ \ (K > 7t) are 
shown. Dashed-dotted lines show eigenfrequencies for a single di- 
rect or complimentary vortex in a LII with one discontinuity. 

K = 7C via the flat phase state. 24 Note that for a < a c the fre- 
quency co + = 1, i.e., belongs to the plasma band. This means 
that strongly coupled vortices behave like a single object with 
a single eigenfrequency co_ . When a — > a c , CO (jc) — ► 0, corre- 
sponding to the instability of the fiat phase solution in favor of 
the AFM statei 20 i 21 i?° For a — > 0, CO- (k) — ► 1 that corresponds 
to the absence of the eigenmodes when discontinuities cancel 
each other and one ends up with a conventional LJJ. In partic- 
ular, for K = K and a < a c (flat phase state) one can calculate 
analytically the eigenvalue as a function of the facet length 
a. Note that the linearized equation describing the eigenvalue 
problem is nothing else but the Schrodinger equation with a 
potential well-- . Some simple algebraic calculations give 

A 2 sin^aVl -A 2 ) + \/l -A 4 cos (a\/l-A 2 ) =0. (16) 

Writing A, = ;C0_, we will obtain co_(7t) as the first root of 
Eq. i )16t . One can show as well that for < a < a c , co_(jt) 
is the only eigenvalue of the system in accordance with the 
reported numerical results above. 



Now we consider a LJJ with two (+K, +k) discontinuities 
with (271 — k, 27t — k) vortices pinned at them, i.e., the \ \ state. 
The same applies to the J J state by substitution K — > 2k — K. 
In simulations we were starting from the two semifluxon state 
|t = \ \ at K = K and were sweeping K towards or towards 
2k. In a symmetric FM state the in-phase mode frequency 
C0 + is smaller than anti -phase mode frequency co_, namely 
co + < coo and co_ > COo- Below we calculate the frequencies 
co± for K c (a) < k < 2k. 

The quantity K^(a) (0 < 1Q < k) is the critical value of 
discontinuity at which the ground state f \ ceases to exist (be- 
comes unstable) and the system switches to another state, e.g., 
| The values of co± (k) calculated numerically are shown 
in Fig. |3] Note that the lower eigenfrequency C0+(k) — > 
when K — ► k£ . 

In the limit of a — > 0, one can get an analytic expression 
of the eigenvalue as a function of K. This is caused by the 
fact that in the limit a — > 0, the two fractional vortices form 
a single vortex with a double topological charge. Hence, in 
this limit the largest eigenfrequency is given by co_ = <Mo(2k) 
or, following the choice of vortices presented in Fig. [3] co_ = 
COo(47C — 2k) with COo is given by Eq. i ll 11 1. In the case of K = 7C, 
C0_ = in agreement with the fact that the two semifluxons 
form an integer fluxon. 

Another property that is known in the limit of a — > is that 
there is no eigenvalue other than co_ (see the calculations pre- 
sented in Sec. ED. The internal mode co + enters the phonon 
band in this limit leaving co+ (tc) = 1 . 
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C. Asymmetric AFM state 

Now we consider two vortices with topological charges — K 
and 271 - K > in the AFM state J \ admitted by (+K, +k)- 
discontinuities. In simulations we were starting from the state 
J,t=J [ at K = 7t and were sweeping K towards or towards 
2k. The calculated frequencies co± for < K < 2tc are shown 
in Fig. |H 

Note that although this state is asymmetric (at least for large 
a) in terms of magnetic flux carried by the vortices, its Cu±(k) 
dependence is symmetric with respect to the K = 7C axis, i.e., 
K>± (k) = co± (k). This state is stable for < K < 2%. 

When the vortices are weakly coupled [a > a c (k) = 
aJ^K), Fig. |Ha)] we indeed have an asymmetric state 24 
and splitting of eigenfrequencies with co_ < coo < C0+. In 
the opposite limit when the coupling is strong [a < a c (k), 
Fig-EJb)], a new collective state is formed^ This state has 
only one eigenfrequency co_, while co + = 1 and belongs to 



semifluxon 
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the plasma band. 

The crossover distance a c ' (k) which separates an asymmet- 
ric state from a collective state is a weak function of K as can 
be seen in Fig.|4] The crossover distance a| ^(K) grows from 
a l c l {n) = 7C/2 « 1.57 to 4 r (0) = a l c l {2%) = 21n(l +y/2) « 
1.7628^ Thus, for 7t/2 < a < 21n(l + \/2) the asymmetric 
flux state exists not in the whole interval of < K < 27C, but 
only in some smaller sub-interval which includes K = 7C. In 
the same time the collective state nucleates around K = and 
K = 2k. At a < k/2 the island of the asymmetric state van- 
ishes and the collective state exists in the whole interval of 
< K < 2tz. Vice versa, at a > 21n(l + \J~T) the asymmetric 
state exists for all < K < 2%. 

Since for K = 7C the symmetric and asymmetric AFM states 
are equivalent to the flat phase state, the same limiting behav- 
ior which was discussed at the end of Sec. MI Al apply [see 
also Eq. dl6> and c.f. Fig. |2). 

It is interesting to follow the smallest eigenfrequency 
G)_(k) at K = K in Fig. In Fig. 0a), decreasing the facet 
length a gives smaller co_ (ft) until the eigenfrequency reaches 
at a = a c that corresponds to the transition to the flat phase 
state. After a passes a c , in the flat phase state, further decrease 
of the facet length increases co_ {%). 

For arbitrary K and a — > the discontinuities (+K, +k) be- 
come a single discontinuity +2k. This case is similar to the 
limiting case we presented in Sec. IIIIBI but the topological 
charge of the fractional vortex is now 2n — 2k. Hence, for 
a — ► the eigenfrequency of the collective state approaches 

co-(k) =coo(27C-2k), (17) 

with tt>o defined by Eq. ( II 11 [see the line marked as 0.0 in 
Fig. lib)]. 

D. Asymmetric FM state 
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FIG. 4: Eigenfrequencies £0_ (k) (solid line) and co+ (k) (dashed line) 
for the asymmetric AFM state J [ calculated numerically for different 
distances a between vortices. In (a) and (b) the coupling between 
the two vortices is strong and weak, respectively. The line marked 
as 0.0 in (b) shows the behavior of the eigenfrequency in the limit 
a — » calculated analytically using Eq. 1171 . Dashed-dotted lines 
show eigenfrequencies for a single direct or complimentary vortex in 
a LJJ with one discontinuity. 



Finally, we consider two vortices +K and 2n — K > in the 
FM state ] \ admitted by two discontinuities — K and +K. In 
simulations we were starting from the state TT=1 t at K = K 
and were sweeping K towards or towards 27C. In this state 
the in-phase mode frequency co + < (Do < (0_. The calculated 
frequencies co± are shown in Fig. [5] Again, the plot is sym- 
metric with respect to K = Jt, i.e., g>± (k) = co^(k). 

The asymmetric FM state has no crossovers and at a — > 
turns itself into a 27C vortex. 

One can also recognize that for a -c 1, (— %,+%)- 
discontinuities in a Josephson junction behave like a micro- 
resistor to an integer 27C-fluxon (see, e.g., Eq. (1.2) of Ref. Bll). 
In the limit a — > 0, the micro-resistor strength e from Ref. 
is equal to 2a. 

IV. CONCLUSIONS 

We have shown that an arbitrary fractional vortex has an 
eigenmode that corresponds to the oscillations of the vortex 
around the phase discontinuity point where it is pinned. We 
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FIG. 5: Eigenfrequencies co_ (k) (solid line) and co+ (k) (dashed line) 
for the asymmetric FM state ] \ calculated numerically for different 
distances a between vortices. Dashed-dotted lines show eigenfre- 
quencies for a single direct or complimentary vortex in a LJJ with 
one discontinuity. 



have derived the eigenfrequency for a direct and a complimen- 
tary vortex as a function of the topological charge (magnetic 
flux) carried by the vortex. The eigenfrequency corresponds 
to a discrete "impurity" energy level within the plasma gap. 

For the case of two coupled vortices we showed that the 
eigenfrequency splits into two modes that correspond to the 
in-phase and anti-phase oscillations of the vortices. Both of 
these eigenfrequencies were calculated numerically for each 
ground state. For a < a c , two AFM ground states form a 
strongly coupled "vortex molecule" which behaves as a sin- 
gle object and has only one eigenmode. 

The obtained results may have enormous impact on the 
study of vortex dynamics in novel Josephson systems, e.g., 
based on SFS junctions or rf-wave superconductors. The 
knowledge of eigenmodes is a key element in designing clas- 
sical or quantum devices based on fractional vortices. In the 
classical domain this may help to avoid parasitic resonance 
phenomena. In the quantum domain, the eigenfrequency de- 
termines the characteristic frequency of, e.g., quantum tunnel- 
ing/flipping processes. 

To detect eigenmodes experimentally one may use spectro- 
scopic methods or resonant excitation. To make spectroscopy 



one should excite the eigenmode by applying an ac bias cur- 
rent or sending an electromagnetic wave along the LJJ from 
the edge. Then one should measure e.g. the critical current as 
a function of the external excitation frequency co. The ac bias 
current can be applied by coupling the junction bias leads to 
a microwave antenna (and avoiding the induction of the mi- 
crowave signal in all other parts, e.g., injectors) or by embed- 
ding the junction into a resonator. If we choose to send the 
microwaves from the edge, the length of the junction should 
not be very large and the frequency should not be very low, so 
that the condition 




< 1 



is satisfied, i.e., the microwave signal should not decay sub- 
stantially while it propagates from the edge to the middle of 
the junction where the vortex (molecule) is created. Spec- 
troscopic methods may help to distinguish between a direct 
K- vortex and its complimentary (2k — k) -vortex by measuring 
the eigenfrequency (except K it). 

Another possibility to detect an eigenmode is to create a 
fractional vortex in an annular LJJ and let the integer fluxon, 
which can be injected if needed^^, run around the LJJ pe- 
riodically colliding with the fractional vortex. If the vor- 
tex's eigenfrequency is a multiple of the collision frequency, 
one should see resonances on the IVC. The first study of a 
fluxon-semifluxon interaction in an annular LJJ was recently 
reported 33 , but the resonances, which appear due to the exci- 
tation of the eigenmodes, were not investigated in detail. 

For future research, it is interesting to study an infinite ar- 
ray of fractional vortices (ID vortex crystal) which has only 
optical branches in the dispersion relation because of the vor- 
tex pinning. The energy bands of such an array can be tuned 
by changing K or bias current y during experiment. 
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Here we consider fractional vortices which appear as a result of 
the phase discontinuity— Other types of arbitrary fractional 
vortices may appear in a LJJ with a strong second harmonic in the 
current-phase relation, which can be present either intrinsically or 
due to a faceted grain-boundary (frequently alternating regions of 
size <C Xj of negative and positive critical current). In contrast to 
such "splintered" vortices— i—i— which are not pinned, the frac- 
tional vortices discussed here are pinned at a K-discontinuity. The 
mixture of both possibilities, i.e., a fractional vortex at the bound- 
ary between a (p-region with the second harmonic and a 0-region 
without the second harmonic is also possibla— 
A fluxon has two trivial eigenvalues, i.e., at zero frequency (with 
multiplicity two) and at the edges of the plasma-band™- The zero 
eigenvalue is due to the translational invariance. The eigenvalue at 
the edges of the plasma-band bifurcates and becomes an internal 
mode in the presence of a perturbative term in the sine-Gordon 
equation. 

The splitting of an eigenvalue out of the edges of the phonon band 
cannot be calculated using the formal expansion method. This is 
an indication of an exponentially small splitting of the eigenvalue. 



